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Abstract 

Employing Forrester-Ha method of Jack polynomials, we derive an integral identity con- 
necting certain A'^-fold coordinate average of the Calogero-Sutherland model with the n-fold 
replica integral. Subsequent analytical continuation to non-integer n leads to asymptotic 
expressions for the (static and dynamic) density-density correlation function of the model as 
well as the Green's function for an arbitrary coupling constant A. 

1 Introduction 

The model of one-dimensional particles interacting via the inverse square potential, known as the 
Calogero-Sutherland model (CSM) has attracted a great deal of attention H, |[ ^ It has 
long been realized [Q that for the three special values of the interaction parameter A = 1/2, 1, 2 
the CSM is closely related to the random matrix theory (RMT) ||] of the three Wigner-Dyson 
ensembles: orthogonal, unitary and symplectic respectively. Recently, using the method of Jack 
polynomials, Forrester Q found the correlation functions of the CSM for any integer A. The 
method was latter generalized by Ha Q to any rational A = p/q, where p and q are mutual primes. 
In the latter case the density-density correlation function of the model may be written as the 
(p + g)-fold integral over momenta of p hole and q particle excitations |j. This seemingly 
makes the answer a highly irregular function of the coupling strength A. 

In this paper we employ the replica method, borrowed from the theory of disordered systems 
and specifically RMT ^, |lOj, to treat the CSM at an arbitrary value of A. The idea is to 
write the particle density operator as the n-th power of the characteristic polynomial in the limit 
n — > 0. In the context of the CSM, this appears to be a particularly simple operator which creates 
at most n holes and no other excitations [^. As a result, for an integer n, the correlation function 
of the characteristic polynomials may be expressed as an n-fold integral over hole momenta for 
any coupling constant, A. 

Mathematically, this statement boils down to the powerful asymptotic (TV — > cxd) identity 
between the iV-fold coordinate average of the n-th power of the characteristic polynomial and 
the n-fold momenta integral. For A = 1 (unitary RMT) such integral representation has been 
obtained earlier by Verbaarschot and Zirnbauer by mapping the problem onto the non-linear 
tr-model. We derive here the identity valid for any A employing the Forrester-Ha method of Jack 
polynomials. As of today we are not aware if this identity has a geometrical or field-theoretical 
derivation apart from A = 1/2, 1, 2. The simplicity of the final expression suggests the possibility 
that such a derivation exists. We stress, however, that for any A the identity has clear quantum- 
mechanical interpretation. We hope also that it may prove to be useful in the theory of the 
Coulomb gases. 

The calculation requires evaluation of the resulting n-fold integral for any integer n. We ac- 
complish this goal in the asymptotic regime a; 3> 1 (a; is the scaled distance) by summing up 
contributions on n + 1 saddle points (n — 1 of these saddle points may be dubbed as "replica 
non-symmetric"). The subsequent analytic continuation, n — > 0, is analogous to the one em- 
ployed recently in the RMT context P, nol . Some mathematical aspects of this procedure are not 
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completely clear at the moment |Ll| . For rational coupling constants A our procedure is in exact 
agreement with the existing results [T^ and allows one to obtain them in a much simpler way. 
In addition it leads to expressions formally valid for any A. 

The paper is organized as follows. In section |^ we define the CSM and introduce the replica 
method. In section || we derive the integral identity discussed above. Section ^ is devoted to the 
calculation of the n-fold integral in the large x limit and analytical continuation, n — > 0. Finally, 
our results and some conclusions are summarized in section |^. Technical details are presented in 
three appendices. 



2 CSM and replica 

2.1 Eigenstates of CSM 

Consider N interacting particles on the unit circle. The CSM is defined by the following 
Hamiltonian: 

H^±i.Af + 2Y^I^^, (1) 

3 = 1 3<i ' ^ 

where Zj — e^^' are the coordinates of the j-th particle, dj = d/dzj and A is the coupling constant. 
The first term is the usual kinetic energy, while the second one describes inverse square interactions. 
The spectrum and the many-body eigenfunctions of the model may be found exactly using the 
theory of the symmetric functions . For the sake of completeness and to set notations we 

present here the basic features of this solution (for more details see Appendix ^) . 
The ground state wave-function of the CSM |l) is jl^ 

(z|0) = *o(z) = [^N{zt , (2) 

where 

1\n{z) = \{{z^~ Zj) (3) 

3<i 

is the N X N Vandermonde determinant and the normalization constant is given by 

A/- = A/-(A,iV) = (2.)-0±^. (4) 

Excited states are labeled by the quantum numbers (essentially momenta) (k) = (ki, K2, . . . , kn)i 
where the non-negative integers Kj are ordered as ki > . . . > kat > 0. In the mathematical liter- 
ature the quantum numbers (k) are called partitions [p^ and represented by the Young diagrams, 
with Kj cells in j-th row. The example of the Young diagram for the partition (k) = (7, 6, 5, 5, 3, 2) 
sometimes written as (7, 6, 5^, 3, 2) is presented in Fig. ^ In our case number of rows in a Young 
diagram cannot exceed N, whereas number of columns is unlimited. The state (k) has well defined 
total momentum, whose value is given by the number of cells in the corresponding Young diagram 
and denoted as |k| = - The coordinate representation wave- function of the excited state 

corresponding to a given partition (k) may be written as 



(z|(«)) = ^^^^iz) = ^)jli){z) X [^^{zt , (5) 

where the symmetric functions J^^^{z) = J^^^{zi, . . . , zjv) are known as Jack polynomials p^ , p!4| . 
Definitions and some properties of Jack polynomials along with the normalization constants, Af(K)j 
are presented in Appendices ^ ^ 
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a) 



b) 



Figure 1: a) Young diagram for partition (k) = (7,6,5,5,3,2) representing an excitation of 
particles in the Calogero-Sutherland model, b) The Young diagram for the dual partition 
(k) — (6, 6, 5, 4, 4, 2, 1) representing the excitation of holes. 



2.2 Correlation functions 

We focus first on the ground state, static density-density correlation function of the CSM: 

Ri0)^Po'm9)p{om, (6) 

where p{9) — ^(6*) — po denotes the fluctuating part of the density operator p{9) — J^j S{6'-6j) and 
Pq — N/2n is the mean density. Due to the translational invariance of the model this correlation 
function depends only on the relative angle between the two points. One may represent the 
fluctuating part of the density operator as the Fourier series: 

where the Fourier coefficients of the density, pk ~ P-k, are given by 

N N 

In the theory of symmetric functions the Fourier coefficients 2T:pk are referred to as power sums 
of variables Zj. Since the Jack polynomials constitute a complete basis in the ring of symmetric 
polynomials the power sums pk can be decomposed into linear combinations of J^^^(z), which is 
the basis of Ha's Q approach. Here we employ a different method. 
Let us define the characteristic polynomial (or vertex) operator as 

N N 

z{0) n (i - ^"'^'"'^0 - n (1 - ^"''-.) • (9) 

In the RMT this object is known as the characteristic polynomial of the unitary matrix with 
eigenvalues Zj or spectral determinant. 
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Rewriting the characteristic polynomial in the form 

N 



Z(0) = exp|^ln(l-e-^%)| =exp|-2^^^^^i^| , (10) 

one notices that the fluctuating part of the density is given by 

p{e) ^ -delm\nZ{0) . (11) 

TT 



To ensure convergence of the series in Eq. ( [LOD one should assume an infinitesimal negative imag- 
inary part of 9: 9^9 — iO. The correlation function Eq. (^) may be expressed therefore through 
the correlation of logarithms of the spectral determinants: 

-2 

R{0) = -^d^Re{0\\nZ(9)lnZ(0m . (12) 
The replica trick is to write the logarithm as 

\nZ= lim - (Z" - 1) . (13) 

n^O n 

As a result one finally obtains for the correlation function 

R{9) = - hm -1^ Re (O|Z"(0)Z"(O)|O) , (14) 

If one was able to calculate the expectation value on the r.h.s. for any real n the limiting procedure 
would be straightforward. The non-trivial aspects of the replica method show up if the calculation 
is possible only for positive pl| integer n. As we shall see in the next section this is indeed the 



3 Correlations of characteristic polynomials 



In this section we calculate the large N limit of the correlation function of the characteristic 
polynomials for integer n. In the coordinate representation it takes the form of the following 
A^-fold integral 



(0|Z"(6l)Z"(0)|0) =AA d^6\/\N{z 



|2A 



N 



e'zj) (1 - z,) 



(15) 



where the normalization constant is given by Eq. (Q). 

Consider the ground state expectation value of the s-th power of the characteristic polynomials, 
where s is an arbitrary real number. Inserting resolution of unity, one finds 



{Q\Z%9)Z^{0m = 5](O|Z^(0)|(«)>((«)|Z^(O)|O) ^ l(WI^'(0)|0)|^e^l'^l' 



(16) 



where the sum runs over all possible partitions (with number of rows less than N). The last 
equality in Eq. (|l^) is a trivial consequence of the translational invariance and the fact that the 
state |(k)) carries momentum Since ^''(0) is a symmetric function of z it may be expanded 
in the basis of Jack polynomials. Employing the orthogonality relation for Jack polynomials, one 
obtains the matrix elements ((k)|^*(0)|0). Details of this calculation are given in Appendix 
The result is 



|Z^(e)Z^(0)|0)=^|A(,)(A,iV)P ({-.} 



(17) 
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where ^(k) (-^j is an s-independent factor, which is given in Appendix H and the symbol {— s}^ 



is defined as 



n (-^-A(j-l) + (a-l)) 



j,ae{K) 



(18) 



Here the product runs over all cells of the Young diagram for a given partition (k) with j = 1, 
being the vertical coordinate and a = 1, 2, . . . the horizontal coordinate of the cell. 

At this point the advantage of the replica approach may be appreciated. One notices that for 
s — n - positive integer, the factor {—n}^^^ defined by Eq. ( p^ ) is zero for any partition (k) which 
has more than n columns in its Young diagram. Indeed, for j = I (the upper row) and a = n -\- 1 
the factor on the r.h.s of Eq. ([l^) nullifies for any A. As a result the set of partitions which 
contribute to the correlation function, Eq. (17), is greatly reduced. It contains only partitions 
which have not more than N rows and not more than n columns in their Young diagrams. 

The similar drastic reduction in the number of contributing partitions takes place in the Ha's 
approach However, in that case it happens only for rational X — p/q and the geometry of 
the relevant Young diagrams is determined by p and q. In the replica approach the fact that n is 
integer takes care of diagrams selection, whereas A may be arbitrary. In other words, the operator 

is a particularly simple one: by acting on the ground state it creates not more than n hole-type 
excitations. As a result, number of partitions (excited states) needed to represent Z"|0) is finite 
and each partition can be parametrized by the heights of its n columns. 

The next step is to take the thermodynamic limit, N oo. To this end we parameterize the 
relevant partitions by the height of the columns in their Young diagrams. This introduces the 
dual partition (k) = (ki, K2, ■ ■ ■ , Hn) whose Young diagram is just a transposition of the diagram 
for the partition (k). Since Young diagram of (k) has no more than N rows and no more than 
n columns, < iia ^ N and a = l,2...n. Physically ka describe momenta of n holes; each 
cannot exceed N since the "Fermi" sea is bound (recall that kj describe unbounded momenta of 
N particles). In the limit N oo one may pass from summations over ka in Eq. ( p^ ) to integration 
over ka = ka/N, where < fca < 1, simultaneously taking the limit in A(K)(A,iV —^ oo). The 
corresponding procedure is described in details in Appendix |c[ As a result one finds in the large 
N hmit: 



(O|Z"(0)Z"(O)|O) ^N-Ci\n) Jd"k\Ar^{k)\ 



where the normalization constant is 



2/A 



Y[ka{l-ka) 



1/A-l 



exp I iN0j2ka \ , (19) 



C(A,n).n^+l/^) 
^ ' ' J-J; T(l + a/ 



11 r(l + a/A)r(a/A) 



(20) 



Finally, employing Eqs. (|T^ and (Q), one may express the result as the identity of the two 
integrals 



lim N-'^Af(X,N) 

N->OD 



27r 



2A 



C(A,n) /d"A:|A„(fc)|2/^ 



N 



n(i 



J2-iTix/N 



Zj) (1 - Zj) 



\{ka{l-ka) 



1/A-l 



exp < ka \ , (21) 



where we have introduced the scaling variable x = po9 = NQ/{2'k) which measures distance in 
units of the mean inter-particle separation, to emphasize that the left hand side has a well defined 
thermodynamic limit. 

For A = 1 (unitary RMT), the identity Eq. (|2^) was derived by Verbaarschot and Zirnbauer 
1^, using fermionic replica pi| and mapping the random matrix problem onto the non-linear 
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(T-model on the coset space U{2n)/U{n) x U{n). In that formulation parameters ka play the role 
of the (compact) angular coordinates on the coset. In a similar way one may derive Eq. ( pl| ) for 
A = 1/2 and 2 (orthogonal and symplectic RMT) It would be extremely interesting to see if 
there is geometrical field-theoretical interpretation of Eq. ( pT| ) for other A. 

The integral identity Eq. ( |2l| ) is exact for any 9. In particular, for a; = its right hand side 
may be evaluated with the help of the Selberg integral Q leading to 

hm 7V-^(0|Z«(0)Z"(0)|0) = n /^"^^^ . : (22) 

air((a + n)/Aj 

result obtained previously by Keating jl^ for A = 1. For x ^ the integral on the r.h.s. can 
not be evaluated exactly. In the next section we shall evaluate it in the asymptotic limit a; 3> 1 
employing the stationary phase method as suggested by Forrester |^ and Yurkevich and Lerner 
0. 



4 Replica Limit 

For X 3> 1 the dominant contribution to the integral on the r.h.s. of Eq. ( ^l| ) comes from the 
end points of the interval. We expand the integrand in I variables about 1 and in remaining n — I 
variables about to leading order. Extending limits of integration to infinity, one obtains for the 
integral: 

J2 f ?) (-2^*^)^; {27Tix)e^^'^' , (23) 



1=0 

where Im are Selberg integrals B 



Q Of — X Q, — X 

As a result, one obtains for the correlation function of the characteristic polynomials in the asymp- 
totic limit N ^ oo, x 1 

{0\Z-{e)Z-{Om=N^Y.[PnW] S E' (25) 

/=0 {~2'!Tix) > {2TTix) > 



where 



I n—l 

n r(l + a/A) n r(l + a/A) i 

KW - (?) „ = (?) n , X (26) 

nr(l + a/A) \^K=ir[l + in^a+l)/X) 



and F^{X) = 1. We have extended summation over / to infinity in Eq. (|25|), since for an integer n 
the binomial coefficient (") is identically zero for ^ > n. For RMT ( A = 1/2, 1, 2 ) Eq. ( p5| ) was 
derived recently from the saddle point treatment of the non-linear cr-model |l^. In that case 
the factor F^{X) is the volume of the degenerate coset space manifold G{n) / G{l)G{n — I), where 
G is an appropriate symmetry group, G = Sp, U,0 for A = 1/2.1.2 correspondingly. 

Up to this moment all our calculations were performed for integer n only. Next we shall 
continue analytically Eq. ( |25| ) to arbitrary real n and then take the replica limit n — s- 0. To this 
end we notice that 



^ ' ^ n-i — ^ , n ^ (27) 



I J r(i + or(i + n-0 
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This expression guarantees cancellation of the factor n ^ in Eq. ( p^ for small n. As a result one 
obtains for I > 1 and n — > 

nm^ntip-(l^_^il0>_,Oi„'). (28) 



Finally, employing Eq. ([IJ) and taking the replica limit, one finds for the correlation function 

^(-) = -^ + (4^7X^(2'^'-)' (29) 
1=1 ^ ' 



where 



n r(l - a/A) 



The first term on the r.h.s. of Eq. (|29| ) comes from the I ~ contribution. It may be called 
replica symmetric, since all n integrals are calculated around ka — (the left "Fermi" point). 



The oscillatory contributions on the r.h.s. of Eq. (29) may be dubbed "replica non-symmetric", 
since in I out of n replica contribution to the integral comes from the vicinity of fca = 1 (the 
right "Fermi" point). 

For rational coupling constant X — p/q the sine function in Eq. (|30| ) vanishes for a = p. As a 
result, di{p/q) — foi I > p and therefore the series in Eq. ( |29| ) contain only p oscillatory com- 
ponents. For this case our results Eqs. (|2^), ( |30| ) coincide exactly with the asymptotic expression 
for the density-density correlation function obtained by Ha [1). In Ref. ||^ the coefficients df{X) 
are given in terms of the Dotsenko-Fateev integrals. One must employ the integral identity of 
Forrester and Zuk |12 to show the equivalence of coefficients in Ref. js) and Eq. (^0|). In fact our 



approach together with Ref. may be considered as an independent evaluation of some of the 
Dotsenko-Fateev integrals. We believe that for irrational A the series in Eq. ( p9| ) is an asymptotic 
one, but its exact meaning and possible re-summation procedure requires more study. 



5 Discussion of the Results 

Let us briefly discuss the generalization of the above procedure to the dynamical correlation func- 
tion R{x,t) = pg ^(0|p(a;, i) /5(0, 0)|0). The characteristic polynomials and repfica are introduced 
in exactly the same way as for the static correlation function. One should take into account that 
operators at different times do not commute and that p-k{t) — p\{i)- As a result one flnds (cf. 
Eqs. (|l|) and (|l|) ) 

i?(x,t) = - hm ^|((K)|Z"(0)|0)|2cos(|^|x/po)e-^*(^"')-^°) , (31) 

where the energies of excited states, are given in Appendix The matrix elements of 

^"(0) and their thermodynamic limit are exactly the same as for the static case. The thermody- 
namic limit of the excitation energy in terms of hole momenta in the proper frame of reference 
is given by Eq. (^). As a result, the integral on the r.h.s. of Eq. (|l9| ) acquires the factor 
exp{— 27ri[Wii ka{l — ka) T ^aW, where Vs = 27rpQA is the sound velocity in the rescaled 
coordinates. 

Asymptotic calculation of the n-fold momenta integral in the limit x ^ 1 and/or Vgt ^ 
1 follows closely the static case. The only difference is that the product of Selberg integrals 
2In-i{—2'jTix)Ii{2iTix) in Eq. (|2^ ) is substituted by /„_;(27rj[wst — x\)Ii{2T:i[vst + x]) -\- {x ^ —x). 
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Performing the analytical continuation n — > and employing Eq. 
dynamic correlation function 



1 



1 



1 



2E- 



one finally finds for the 



cos{2tiIx) , (32) 



where the coefficients c?;(A) are given by Eq. (|30|). 

The method also enables one to calculate the single particle Green function of the model. As 
it was shown in Ref. the annihilation operator of particle acting on the ground state of + 1 
particles is proportional to the action of on the N particle ground state, i.e. 



4-^(0;, t)*(0, 0)|0) = (A + |((«;)|Z^(0)|0)pe'(2-l«l/^— 

J\l (A, iV j — ^ 



(33) 



Therefore to evaluate the Green function one has to perform analytical continuation of 
integer n to n = A. This yields in the large A^ limit 



5f) from 



(0|*+(x,t)*(0,0)|0) 



1 



1 



271'^'-"^^"^-' \2'ni{vst — x) J \2Tri{Vst- 

(34) 

where the coefficients Fj,(A) are still given by ( |2^ ) with n A. For rational X — p/q they vanish 
for / > p as in the case of the density-density correlations, so the sum in (^) contains a finite 
number of terms. 

In summary, the correlation functions of the integer powers n of the spectral determinants may 
be evaluated as a sum over partitions with no more than n columns. In the thermodynamic limit 
such a sum takes a form of the n-fold integral, which may be evaluated for the large inter-particle 
separation. The subsequent analytical continuation to non-integer n results in the density-density 
correlation function (n — s- 0) or in the single particle Green function (n —f A). These calculations 
do not require rational values of the coupling constant A and are formally valid for any A. The 
final expressions, however, take the form of divergent series for irrational A, whose exact meaning 
is not clear at the moment. For rational A the results coincide with those of Ref. 1^. 

Discussions with B. Altshuler initiated this project. We are also indebted to Y. Avron, J. 
Feinberg, S. Fishman, B. Simons M. Zirnbauer and A. Ludwig for valuable conversations. This 
research was supported in part by the U.S. -Israel Binational Science Foundation (BSF), by the 
Minerva Center for Non-linear Physics of Complex Systems, by the Israel Science Foundation, by 
the Niedersachsen Ministry of Science (Germany) and by the Fund for Promotion of Research at 
the Technion. 



A Eigenstates of the Calogero-Sutherland model 

Consider action of the kinetic term of the Hamiltonian Eq. (^ on the ground state, Eq. 
is convenient to rewrite as exp{Aln A(z)} : 



N N r 

J2 i^Af [A(^)]' = E A {z,d,f In A(z) + X'[z,d, In A(z))^ 



The logarithmic derivatives are evaluated as follows: 



ZjdjliiA{z) = ^ 



{ZjdjflllAiz) 



[A(z)]^ 



), which 



(35) 



(36) 
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As a result the action of kinetic term Eq. (|3^) is equivalent to the multiphcation by the following 
factor 

A(A-l)Er^ + ^^f E 7 i , + (37) 

^^(z.-z,)^ l^^^^,(z.-z,)(z.-z.) 2 j 

The first term in this expression cancels exactly the action of the interaction term of the Hamilto- 
nian Eq. ([|), while the second term may be easily shown to be a constant. That proves that the 
wave-function Eq. (||) is indeed an eigenstate of the Hamiltonian Eq. (|l|) with the energy 

y 3 2 y 6 

which may be considered as a sum of N independent one-particle energies ej — }?{N ~ jY 
for j running from 1 to N . The normalization of the ground state wave function is obtained 
straightforwardly with the help of the Selberg integral j|] 

"^""^ ni-^' ,.2A_ r(l-HAiV) 
(2^)ivlil^ ^ I -r^(l + A) ^^^^ 

resulting in the expression (||) for the normalized ground state. 

To construct the excited states of the CSM one multiplies the ground state wave function by 
some symmetric polynomials ^{z) J[z)"^q{z). The action of the kinetic term on this wave 
fimction generates the term (^7|) through the action of derivatives on ln^'o(2:) and the term, which 
involves the derivatives of In J(z), leading to the eigenvalue equation: 

(i/o + Ai/i) J{z) ^(e~ Eo) J{z) , (40) 

where 

Ho = J2{zA)\ i7i=E^^(z,9,-z,a,). (41) 

. , ■ Zi Zj 

Solutions to this differential equation J(z) = J^^^{z) are the symmetric Jack polynomials 
labeled by the quantum numbers (k) — (ki, K2, ■ ■ ■ kn)- These numbers are non-negative integers 
organized in the non- increasing sequence Ki>K2>...>t7v>0 are known as parts of the 
partition (k). Since Jack polynomials are homogeneous functions, under translations Zj = e*^^ 
uzi — e*(^+*j) they transform as 

J'''{uz)^u\^\j'''{z). (42) 



(k) V ^ (k) 

The total momentum of an excited state is therefore 



Pi^t^Po + l^l, (43) 

where Pq = XN{N — l)/2 is the total momentum of the ground state and |k| = ki + K2 + ■ ■ ■ + kn- 
The energy of the excitation (k) is given by 



AT 



JV 

^0 + E ('^^ + 2Ak, (N - j)) . (44) 



Each term in this sum can be represented as an energy gain ^A(A^ — j) + Hj^ — (^X{N — j)j of 

particle excited by momentum nj from its ground state momentum A(A^ — j). 

It is customary to represent partitions (k) by Young diagrams, a graph with Hj cells in the 
j-th row. A dual partition is obtained by using the numbers of cells, ka, in the a-th column of the 
Young diagram. 
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Figure 2: Another graphical representation of an excited state depicted in the previous figure. The 
excitations of particles and holes from their ground state position are related through the duality 
of the partitions. The number of particles is iV = 8 and the coupling constant A = 3. 



For example, the partition dual to (k) = (7,6,5,5,3,2) ( the one depicted in Fig. |l] ) is 
(k) = (6,6,5,4,4,2,1). It can be represented by transposing the Young diagram for (k). The 
dual partition can be given a nice interpretation in terms of hole excitations. Consider another 
graphical representation of the eigenstates of the Calogero-Sutherland model shown in Fig. 
Each particle in the ground state is represented by black box of length A, which corresponds to 
the momentum X{N ~ j) added to the ground state. The excited state described by partition (k) 
is represented by moving j-th particle (j = 1 corresponds to the uppermost particle) by Kj units 
up. Some holes of integer length will appear between the particles. Recording the shifts of the 
holes from their ground state positions, one finds that they are described by the partition (k) dual 
to (k) as depicted in Fig. |^. In terms of the dual partition [k] the energy ( |4^ ) of the excitation 
reads 

^(«) " = ^ (^^'*" + " + - ^^a^ ' (45) 

while the total momentum is given by P(„) ^ Pq — \k\, since the number of cells can be counted 
row- wise as well as column-wise. In the limit of large N the leading behavior of Eq. ( ^5|) is obtained 
by rescaling ka — Nka and reads — Eq N"^ -^^0(2 — ka)- Here we must recall that we 
have chosen to work with the ground state with non-zero total momentum. That is equivalent 
to describing a system with zero-momentum ground state from the moving frame. Translation 
back to the frame where the ground state has zero momentum changes (only) expression for the 
excitation energy to 

S(K)-^0^iV2^Afca(l-fca). (46) 
a 

This is the expression we use to calculate the dynamical correlation function, which belongs 
therefore to the rest frame of the ground state. 
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The correspondence between dual partition and hole excitations is a non trivial fact and is 
crucial in our approach to the Calogero-Sutherland model. 



B Matrix Elements 

Here we describe how various scalar products between symmetric functions can be represented in 
terms of partitions. Let us denote hy j = 1, . . . , N the vertical coordinate and by a = 1, 2, . . . the 
horizontal coordinate of the cell in a Young diagram. The cell (1,1) is the upper left corner of the 
partition. Then the symbols { }^^^ and [ ]^^^ are defined by the following product over all cells in 
the Young diagram for (k): 

= n (^-A(,7-l) + (a-l)) (47) 
\s]l, ^ n (5 + («-l)A-(j-l)). (48) 



The Jack polynomials are orthogonal with respect to the scalar product 

TTI iOi ie,|2A 1/A._. .i/A. . ^ 



where the factor j^^ may be written as the product over the cells of the Young diagram: 



•?W= n (^a"J + («:j-a)/A + l/A)(«a-J + («, -a)/A + l) . (50) 

j,ae{K,) 

Employing these expressions, one finds for the normalization constant J^(k) of excited states of the 
CSM (cf. Eq. (|)): 

Jack polynomials form an orthogonal basis in the vector space of symmetric functions and can 
be used to expand any symmetric function. For example, the power sums, representing the Fourier 
coefficients of the density operator, Eq. (||), are given by the following linear combination of Jack 
polynomials: 



k ^ [0' 



/lA 

^1/A, 



E4 = iE-F^<)'(-)- (52) 

j \K\=k -Hk) 

The prime in [0']^^^ denotes that the cell (1,1) is not included in the product in Eq. (^8|). More 
important for our discussion is the expansion of any power of the characteristic polynomials: 

z%o) = Yin - z,y = ^-^J'd)^^) • (53) 

3 (K) ^ -^W 

In terms of normalized wave functions of the CSM, Eq. (|^), the action of the s-power of the 
characteristic polynomial operator on the ground state may be written as the following linear 
combination: 

Z''(f?)|0) =E%)(^'^){-4(«)e~^'''"l('^))' (54) 
(«) 

where 
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C The large A'^ limit 

In this appendix we evaluate the thermodynamic hmit N oo for the correlation function of 
the characteristic polynomials, Eq. (17). To this end we need to consider the partition-dependent 
coefficient 



(56) 



Let us start with {— "}('^) given by the product (^^ over cells of the Young diagram representing 
(k). It is clear from this definition that {—n}'^^^ = if ki > n. One can choose the heights 
0<'^a<-^; a=l,...,nof the first n columns as the coordinates, since only the partitions with 
at most n columns contribute. In these coordinates one can rewrite 



(.) - (-A)I'^I n n ((^^ - « + 1)/^ + J' - l) = i-^t^ J{[{n-a + 1)/A 

a— Ij— 1 a— 1 

Here we have used the definition of the Pochhammer's symbol: 

T{z + n) 



(z)„ EE z{z + l)(z + 2) . . . (z + n - 1) = 
As a result, one finds 



A2|k| ]J 



r(z) 

V'(ka + {n-a+l)/\^ 
r2((n-a + l)/A) 



(57) 



(58) 



(59) 



The factors and [A^ + 1/A — 1]^^-) are calculated in the same way by writing the products, 

Eq. (E^), column- wise: 



mi 



a=lj=l 

n ka 

[7v+i/A-i]f,)=nn(^+«/A-j 



a=lj=l 



r(^iV + a/A+ 1 - 1/a) 
5ir(A^-K, + a/A + l-l/A) 
" r(7V + a/A) 
a=ir(^N-ka + a/X^ ' 



(60) 



(61) 



Next we calculate the factor j^^^ defined in Eq. (|5^). We divide the Young diagram representing 
the partition (k) into n rectangles having the fixed width b — 1,2, ... ,n and height kh — Ht+i as 
shown in Fig. g. Note that when kb = the rectangle b has zero height. With these notations 
the expression (EO) can be written as a product over the rectangles: 



nn n («a-i + (6-a + l)/A)(«,-j + l + (6-a)/A 

(jia - Kfc+1 + (5 - a + 1)/A^ r^Ka - Kfe+i + 1 + (6 - a)/A^ 



n b Y 

6=1 a=l r(^Ka - K6 + (6 - a + 1)/A^ T (^ka - Kb + 1 + {b - a) / Xj 



(62) 



Renaming the running index 6 + 1 — > 6 and changing the order of the products, the numerator of 
the last expression may be written as 

n 

Y[ ^(jia - fin+1 + (n + 1 - a)/A^ ^(jia - Hn+i + 1 + (n - a)/A 



a=l 
n 

X Y[ r(ka-fib-{a-b)/X^r(ka~iib-{a~b)/\ + l-l/\ 

b=a+l 



(63) 
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Figure 3: A Young diagram for n = 12 divided in rectangles of fixed width Rectangles of zero 
height / = 3, 6, 7, 11, 12 are represented by dots. 
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In the denominator of Eq. ( p^ ) the order of the products is interchanged and the terms with a ^ b 
are factored out, the resuh is 

n n 

r"(l/A) n Yl r(^ka-kb-{a- b)/X + 1/a) r(^ka -kb-ia~ b)/X + l) . (64) 

a=l b=a+l 

Finally, using the fact that k„+i — 0, one obtains 



1 " 

f;r(I7A) n/('^" + + r(^, + 1 + (n - a)/X 

» rf^a - Kb - (a - b)/x) rfka -kb-ia- b)/X + 1 - 1/a) 

n -7 \ f (■ 

b=a+i r( Ka - Kb - (a - 6)/A + 1/A 1 Tl - Kft - (a - 6)/A + Ij 



Combining the above expressions together one finally obtains for the coefficient Eq. 

r(i/A) 



n 

a=l 
n 

a=l 



r2((n-a+l)/A) 

r(Ka + (?i-a)/A + l/A) t{n -ka + a/x}j t{n + a/X + l - l/x}j 



r 



(Ka + (n-a)/A + l) r(A^-Ka + a/A + l-l/A) r(iV + a/A) 



- rf^, - Kb - (a - 6)/A + 1/a) r(K, - Kb - (a - 6)/A + l) 

IT — ^ ■ (66) 

r(K, - Kb - (a - 6)/a) r^Ka - Kb - (a - 6)/A + 1 - 1/a) 



X 

b=a+l 



The thermodynamic limit, — > oo, is obtained with the help of the following asymptotic relation: 

|z|^oo 1 (^) 

Rescaling the column heights as ka — Nka, one obtains for the thermodynamic limit of the 
expression (pq) 



. ^^'^^^ — r n ik^-hr^'fiki/^-\i-k.)'/'-'. (68) 

a=l (n - a + 1)/Aj i<a<b<n a=l 



TV 



We relax the requirement /ci > fc2 > . . . > fc„ allowing all the variables to run independently 
between and 1. To do so we have to divide by n\ and take the absolute value of the Vandermonde 
determinant: 

n (fca-fcb)2/^-> |A„(fc)|2/\ (69) 

l<a<b<n 

Finally, changing the sum over partitions to the integrals over ka and absorbing the factor A^^" 
into the integration measure, one recovers Eqs. ( |l9| ) and (pO|). 
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